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A fundamental limit to the sensitivity of optical interferometers is imposed by Brownian thermal 
fluctuations of the mirrors' surfaces. This thermal noise can be reduced by using larger beams which 
"average out" the random fluctuations of the surfaces. It has been proposed previously that wider, 
higher-order Laguerre-Gaussian modes can be used to exploit this effect. In this article, we show 
that susceptibility to spatial imperfections of the mirrors' surfaces limits the effectiveness of this 
approach in interferometers used for gravitational-wave detection. Possible methods of reducing this 
susceptibility are also discussed. 
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I. INTRODUCTION 

Long-baseline laser-interferometer gravitational-wave 
detectors, such as those used in LIGO [1], VIRGO [2], 
GEO600 [3] and LCGT [4], use Michelson interferome- 
try to measure tiny differential changes in arm length 
induced by gravitational waves. Spurious motions of a 
mirror's surface, such as those caused by seismic, ther- 
mal, and radiation-pressure fluctuations, can compromise 
the sensitivity to gravitational wave signals. Brownian 
thermal noise in the dielectric mirror coatings, or coat- 
ing Brownian noise, is known to be the dominant noise 
source in the intermediate frequency band of Advanced 
LIGO [5] and other similar interferometers. 

As described by the Fluctuation-Dissipation Theo- 
rem [6, 7], dissipation via internal friction in the dielectric 
coatings must lead to fluctuations in the thickness of the 
coatings. When the beam spot size is much larger than 
the coating thickness, coating Brownian noise at differ- 
ent locations on the mirror's surface can be considered 
to be uncorrelated. This leads to the following scaling 
law [8, 9], 




which describes how the power spectrum of observed 
coating Brownian noise Sx depends on the intensity pro- 
file I{r) of the optical field which is used to read out 
the mirror motion; i.e. the coating Brownian noise power 
spectrum is inversely proportional to the effective area of 
the optical mode. 

Three families of optical modes have so far been con- 
sidered for mitigating coating thermal noise (see Table I) . 
Among these modes, only the higher-order Laguerre- 
Gauss mode, LGa^s, can be supported by optical cavi- 
ties employing standard spherical mirrors. Due to the 
practical advantages associated with the use of spherical 
mirrors, experimental testing of LGs^a modes has begun. 



Table I: Beam shapes that have been considered for use 
in gravitational- wave detectors, mirror shapes that support 
them and their thermal-noise suppression factors (in power) 
for Advanced LIGO parameters (cavity length L = 4 km, 
mirror radius of 17 cm). 



Mode 


Mirror Shape 


Suppression Factor 


Ref. 


LG3,3 


Spherical 


L61 


[10] 


Mesa 


Sombrero 


L53 


[11, 12] 


Conical 


Conical 


2.30 


[13] 



It has thus far been demonstrated that these modes can 
be generated with high efficiency and resonated in table- 
top cavities with small mirrors [14, 15]. 

An unpleasant property of higher-order LG modes is 
that each LGp^i mode is 2p -\- \l\ + 1-fold degenerate, 
the LGs^s mode being 10-fold degenerate. Mirror fig- 
ure errors will inevitably split each formerly degenerate 
mode into 2p+ |/| + 1 single modes with eigenfrequencies 
which depend on the particulars of the figure error. By 
contrast, a non-degenerate mode, under the same figure 
error, will usually remain as a single, weakly-perturbed 
non-degenerate mode. 

In this work we explore the effects of LGs^a modal de- 
generacy quantitatively, via both numerical and analyt- 
ical methods. Guided by experience with existing inter- 
ferometers we have selected contrast defect as our metric 
of interferometer performance. 

To ground our investigation in reality we incorporate 
mirror figure errors derived from measurements of the 
first Advanced LIGO optics. The creation of these maps 
is described in Section II. 

In Section III, we use perturbation theory to analyze 
the effect of such mirror perturbations on the degener- 
ate subspace which includes the LGa^s mode. We then 
use this newly perturbed set of modes to calculate the 
contrast degradation of a single Fabry-Perot arm cavity 
analytically in Section IV A. 

In Section IV B, we utilize a sophisticated numerical 
field propagation code to confirm analytical results and 



2 



examine a more complicated interferometer topology. 

In Sec. V, we explore two methods of mitigating con- 
trast degradation; neither of the methods was ultimately 
successful. 



II. 
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Figure 1: Power spectral densities of uncoated mirror surface 
roughness. The dashed lines are the measured spectra of three 
Advanced LIGO arm cavity mirrors. A model approximating 
these spectra (black trace, see (2)) was created to generate 
the random maps used in our work. 



In this work we investigate the consequences of realistic 
mirror imperfections on the performance of the LGa^s 
mode. The parameters of the particular imperfections 
applied are therefore significant. 

Fig. 1 illustrates measured surface roughness power 
spectral densities of selected Advanced LIGO mirror sub- 
strates, prior to the application of dielectric mirror coat- 
ings (the PSD plots end up to spatial wavelength of 2 
mm). Based on the measured Initial LIGO optics and 
small optics of Advanced LIGO, we construct an analyt- 
ical model (solid black line) which falls roughly in the 
middle of Advanced LIGO test mass PSDs, 



S{f) cx (1 + (0.04/)' 



(2) 



This one-dimensional function was used to generate ran- 
dom mirror maps which are statistically similar to those 
one might find in an advanced gravitational-wave inter- 
ferometer. Such random mirror maps, used in all as- 
pects of this investigation, were constructed by multiply- 
ing each point of the amplitude spectral density's mag- 
nitude by a random complex number a^ih before trans- 
forming back to coordinate space and appropriately scal- 
ing the result to yield the desired RMS. Scalars a and h 



are drawn independently from a normal distribution with 
zero mean and a standard deviation of one [16]. 

The entire surface is fit by Zernike polynomials and the 
terms corresponding to Piston, tilt and power (Zernike 
polynomials Zg, Z^^ Z^) were removed from our maps 
(the Piston term is irrelevant because the lock process 
adjusts the microscopic length; the tilt term is removed 
to represent the alignment control; the ROC of the gen- 
erated surface is corrected by hand). The RMS values 
quoted are calculated after this subtraction. 

Fig. 2 shows the surface figure of one map generated 
using our algorithm. This map is typical of a larger pop- 
ulation and was selected as a reference to be used in all 
analytical calculations. 
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Figure 2: Surface figure (in nm) of a typical generated phase 
map with Piston, tilt, power and astigmatism terms sub- 
tracted. 



III. DEGENERATE PERTURBATION-THEORY 
ANALYSIS 

A. Laugerre- Gauss modes 

The Laguerre-Gauss modes (LGp^i) are a set of circu- 
larly symmetric modes which can be written in cylindri- 
cal coordinates as [17] 



Up^i{r,(j),z) 



2p\ 
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where uj{z) is the beam radius, tlj{z) is the Gouy phase, 
and R{z) is phase front curvature of the beam. Lp^{x) 
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Figure 3: Normalized intensity distribution of the LG3,3 mode 
at the mirror position (cjq = 0.021 m, 2; = 1997.25 m). 



B. Application of Degenerate Perturbation Theory 
to the Perturbed Fabry-Perot Cavity 

The combination of eigenmodes excited in a cavity de- 
pends on the composition of the incident field and on the 
properties of the cavity itself. In this section we discuss 
how first-order perturbation theory can be applied to this 
problem. We first explore how mirror figure error breaks 
the degeneracy of LG cavity modes before describing the 
phase shift each mode experiences in an optical cavity 
and finally constructing the total field (prompt plus leak- 
age) reflected from a perturbed Fabry-Perot resonator. 



1. Mode Splitting 

Fig. 4 illustrates light propagation in a simple Fabry- 
Perot cavity, introducing the notation employed in our 
formalism. 



Table II: Suppression factors of thermal noise (in power spec- 
tral density) for LG modes with a fixed clipping loss of 1 ppm. 
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Figure 4: Fabry-Perot cavity with a perturbation 6h{x,y) on 
the end mirror. 



is the associated Laguerre polynomial where p>0 and / 
are the radial and azimuthal indices respectively. 

The mode selectivity of the cavity is determined by 
the cavity finesse and the mode dependent phase shift 
(2p+ |/| + l)tlj{z). From this we see that the LGp i mode 
has 2p -\- \l\ + 1 degenerate eigenmodes. For example, 
LGa^a belongs to a 10-fold degenerate space, which can 
be spanned by: LGs^^a, LGo,±9, LGi,±7, LG2,±5 and 
LG4,±i. 

Coating Brownian noise power is proportional to the 
integral of beam intensity, as the scaling law (Eq. (1)) in- 
dicates. In Table II, we present theoretical thermal noise 
suppression factors for selected LG modes. To permit a 
fair comparison, the widths of all modes considered here 
and henceforth were chosen to be 0.018 m, which yield a 
clipping loss [18], due to the finite size of the cavity mir- 
rors, of around 1 ppm. We see that the LGa^a, considered 
by many as the leading candidate for use in gravitational 
wave interferometers, offers a theoretical thermal noise 
reduction factor of ~1.6 compared to a standard Gaus- 
sian beam (LGo,o)- The transverse intensity distribution 
of the LGa^s mode is presented in Fig. 3. 



We use the following standard method for propagating 
cavity fields, 

^r{f) = J K{r,r')^i{r')d^r\ (4) 

or 1^^) = i^l^/), where and are the electric fields 
near the right and left mirror respectively and K is the 
field propagator from the left mirror to the right 

mn = ^e-^i^^-^^'r (5) 

Therefore, if is an eigenmode of the cavity, 

IV^) =RkR'k\il;), (6) 

where R is the reflection operator of the left mirror and 
R! is for the right. Hence, |?/^) is an eigenmode of the 
operator {RKR'K). 

We assume that the mirror on the left is ideal and 
study the consequences of applying a surface figure per- 
turbation Sh to the right hand (end mirror) optic. The 
reflection operator R' can then be written as 

R = Re^'^^^ ^ R{1 + 2ikSh). (7) 
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To obtain the real cavity eigenmodes, we need to solve for 
the eigenfunctions of [RKR'K). Note that the original 
LG33 mode is an eigenfunction of an unperturbed cavity: 



|33) = {RkRk)\33). 



(8) 



As introduced above, the LG33 mode is degenerate 
with 9 other modes, each having eigenfrequency ujq^ thus 
RKRK also has 10 degenerate eigenmodes at uoq. 

For reasonable parameters, modes outside of this de- 
generate sub-space are far enough from resonance that we 
may ignore them in this first-order analysis. Therefore, 
to good approximation, we can assume the new eigen- 
modes of the perturbed cavity are still members of the 
Hilbert space of the original 10 LG modes, which we rep- 
resent by (i=l,2,. . . ,10). These new eigenmodes, de- 
noted are the eigenvectors of the matrix with ele- 
ments {i\RkShRk\j) = {i\Sh\j), {ij = 1,2,..., 10) 



/ {l\Sh\l) ... {l\5h\10) 



(9) 



V(10|(^/l|l) ... (10|^/l|10)y 
Denoting the electric field of \i) as we write 

{i\Sh\j) = J j (t)1{x,y)6h{x,y)(t)j{x,y)dxdy. (10) 

The frequency shift of the degenerate modes introduced 
by the perturbation is then proportional to the eigenval- 
ues of the matrix 



L0i=^l^'\5h\i'), 



(11) 



where k = 27r/A is the optical wavenumber, c is the speed 
of light and L is the cavity length. 
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Figure 5: Frequency shift of LG modes introduced as a result 
of realistic mirror perturbations. 



This quantity was evaluated using the reference map 
shown above (Fig. 2). The RMS roughness of the refer- 
ence was scaled to be similar to those of measured Ad- 
vanced LIGO mirror surfaces (0.3 nm RMS). Results are 
presented in Fig. 5. Here the frequency splits are given 
in Hertz (cj/27r). The frequency shifts are one order of 
magnitude smaller than the advanced gravitational- wave 
interferometer's cavity linewidth. Thus, multiple per- 
turbed eigenmodes will be partially resonant, radically 
distorting the shape of the output field. 



2. The Modal Input- Output Equation 

Above we have shown that mirror figure errors will 
lift the modal degeneracy and split the degenerate LG3^3 
space into distinct states with unique eigenfrequencies. 
We now consider how each of these modes interacts with 
a cavity. 

For any mode |in) injected into an ideal cavity, there 
exists a frequency dependent phase shift between the in- 
put and total refiected or output fields. This can be writ- 
ten as [19] 



7c + 2(a;-a;o)|. \ 
out = r m , 

7c - - ^0) 



(12) 



where 00 is the frequency of the injected field, loq is 
the resonant frequency of the cavity closest to uj and 
7c = cTinput/4i^ is the cavity pole frequency in Hertz. 
Here, Tinput denotes the power transmissivity of the cav- 
ity input mirror; the transmissivity of the end mirror is 
assumed to be zero. 

Suppose that this ideal cavity hosts an A/'-fold degen- 
erate space and that we inject an input mode |in) which 
belongs to this space. If the A^-fold degeneracy is broken 
by some mirror figure error, the new eigenmodes can be 
approximated by |n), n = 1, 2, . . . , A^, where each mode 
still belongs to the original sub- space, but has a new 
eigenfrequency Un- As we shall see in Appendix A, this 
is justified as long as the cavity finesse is high enough 
and the eigenfrequencies of the non-degenerate modes are 
well-separated from this sub-space. 

The output from the perturbed cavity can be obtained 
by projecting the input mode |in) onto the new basis 
\n) and calculating the phase shifts using the following 
relation: 



N 



I out) 



7c + ^(^-^nO| /\/ /I- \ 



(13) 



This procedure was applied to the LG3,3 degener- 
ate space using the resonant frequencies of the previ- 
ous section, Un- With the reference phase map scaled 
to 0.3 nm RMS, the intensity profile of the resulting out- 
put mode is shown in Fig. 6. 
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Figure 6: Intensity distributions of the total field reflected 
from an arm cavity whose end mirror was perturbed by the 
reference map. These distributions were calculated indepen- 
dently via two different techniques. Top - Analytic method 
described in Section IIIB. Bottom - FFT-based numerical 
simulation (see Section IV B). 



IV. CONTRAST DEFECT 

In gravitational- wave interferometers, contrast defect 
is defined as the ratio of the minimum possible optical 
power at the anti-symmetric (dark) port to the power 
incident on the beamsplitter. This quantity can be ex- 
pressed as 



C 



Pas 



(14) 



where X and Y 



Pn = J (^g^t)*^g^t dxdy and ^^s = ^^^^ - 



are labels for the two arm cavities, 
.X -^Y i-ep- 

resents the field at the anti-symmetric port (AS) of the 
interferometer. 

In principle, the dark port could be completely dark. 
However, the presence of intentional imbalances in the 
arms (finite beamsplitter size, Schnupp asymmetry, etc.) 
and unintentional imperfections (mirror shape, scatter 
loss, mirror motion, etc.) result in imperfect destructive 



interference between the fields which recombine at the 
beamsplitter. This imperfect interference leads to the 
leakage of some 'junk' light to the dark port where the 
gravitational- wave signal is also detected. Excess light at 
the dark port can lead to a degradation of sensitivity via 
several mechanisms and compromise the robust opera- 
tion of interferometer longitudinal and alignment control 
systems [20, 21]. Contrast defect is thus a useful metric 
to employ when comparing interferometer configurations. 

The above perturbation analysis shows that the fields 
resonating in the arm cavities of a real interferometer will 
no longer be pure LGs^a modes. Further, the relative 
amplitudes of the quasi-degenerate modes are strongly 
dependent on mirror properties, which will, in general, 
be different for each arm. Hence the perturbed arm cav- 
ity fields will interfere imperfectly at the beamsplitter. 
We therefore expect an LG3 3 interferometer to exhibit a 
larger contrast defect than e.g. an LGo,o mode. We now 
test this hypothesis by analytical and numerical means. 



A. Analytic calculation 

According to Appendix A, the contrast defect in an 
interferometer with one perfect arm cavity and one per- 
turbed arm cavity, can be analytically written as 



(in|out)|. 



(15) 



when the mirror perturbations are sufficiently small 

Appendix A also shows that when a frequency shift 
is small compared with the line width of the cavity, i.e. 
uj — ujn <C 7, Eq. 15 can be approximated as 



e = 



^{33\Sh\n'){n'\Sh\33) - {33\5h\33)^ 



AT 



where 



{33\Sh\pl) = J Sh{ul^Upi) dx dy. 



Alternatively, this can be written as 

2 



€ = 



(8\/27r) 
AT 



where 



(16) 
(17) 

(18) 
(19) 



2p+U| = 10 
(p, 07^(3, 3) 



This suggests that in order to minimize contrast defect, 
one must strive to suppress the projection of Sh onto 
the 9 complex basis functions, given by u^^Upi. Since 
these functions are complex and Sh is real, this actually 
corresponds to 18 basis functions. 
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Figure 7: Contrast defect with the ETM of one cavity per- 
turbed by rescaling the reference phase map: Solid red is from 
the analytical calculation, blue marker is from the FFT cal- 
culation. 



Based on Eq. (13), the contrast defect has been eval- 
uated analytically in the case that only the end mir- 
ror (ETM) of one cavity is perturbed with the reference 
phase map of Fig .2. The perturbation for different RMS 
values was obtained by rescaling the reference phase map. 
The results are shown in Fig. 7. 

One can show that, for small perturbations, perturb- 
ing two cavity mirrors using phase maps derived from 
the same power spectral density function, will, on aver- 
age, result in twice the contrast defect when compared to 
perturbations of a single mirror. However, depending on 
the spatial correlations between the mirrors, the contrast 
defect can be as much as twice the average in some cases. 



B. Numerical calculation 

To confirm the results obtained via perturbation the- 
ory and to extend our analysis to more complicated con- 
figurations, a parallel investigation was carried out us- 
ing numerical methods. We utilized an FFT-based field 
propagation tool - the Stationary Interferometer Simu- 
lation (SIS) [23]. 

SIS is predominantly used to inform the design of the 
Advanced LIGO interferometers and is under continuous 
development at Caltech's LIGO Laboratory. SIS employs 
an iterative procedure to find the stationary fields for a 
given optical configuration and input beam. Mirror sur- 
face maps can be generated from user-defined power spec- 
tral density functions, allowing one to study the effects of 
various hypothetical mirror aberrations. Cavity systems 
are 'locked' using a Pound-Drever-Hall signal [24] to re- 
alize an operating condition similar to that which would 
be observed experimentally. 



Figure 8: Ratio, numerical/analytical, of single-cavity con- 
trast defects calculated with the end mirror perturbed by a 
0.3 nm RMS figure error. 



SIS was used to model Advanced LIGO Fabry- 
Perot arm cavities supporting LGo,0 7 LGa^s and nearly- 
concentric mesa modes [11, 12]. The parameters of the 
LGs^s and mesa cavities (see Table III) were adjusted to 
yield systems with round-trip diffraction loss equivalent 
to that of the fiducial LGo,o resonator (< 1 ppm). In each 
case the input beam remained fixed as the beam which 
was ideally coupled to an unperturbed cavity. To sup- 
press the aliasing effect, a larger FFT grid of 1024x1024 
points ona0.7mx0.7m square was used for all modes. 

Table III: Cavity parameters used in the numerical simula- 
tions. All three resonators had a length of 3994.5 m. The 
mesa radii of curvature refer to the fiducial sphere from which 
the mesa correction profile is subtracted (c.f. [11]). Cavity g 
factor is not well defined for mesa modes. 





LGo,o 


LG3,3 


Mesa 


Ritm 


1934 m 


2857 m 


1997.25 m 


Retm 


2245 m 


2857 m 


1997.25 m 


Cavity g factor 


0.83 


0.16 





Initially SIS was used to simulate a configuration iden- 
tical to that studied analytically. Under these conditions 
both sets of results are in good agreement (better than 
10%) (see Figs. 6, 7 and 8). The flexibility of SIS was then 
used to consider more complex simulations where both 
cavity mirrors are perturbed and to emulate a Fabry- 
Perot Michelson interferometer. SIS was also employed 
to study high-RMS cases in which the analytic approxi- 
mation breaks down. 

For each simulation run, two random surfaces, with 
a specified RMS roughness and a spatial spectrum ap- 
proximating that of the first Advanced LIGO mirrors 
(see Fig. 1), were generated and added to the profiles 



7 



of the cavity mirrors. SIS then evaluated the field 
reflected from the cavity at its operating point, which 
is chosen from the cavity which is locked by using the 
Pound-Drever-Hall error signal. 

Results from two discrete, single-cavity simulations, 
representing the X and Y arms of an interferometer, were 
then combined according to Eq. 14 to estimate interfer- 
ometer contrast defect. Multiple trials were conducted 
at each value of RMS surface roughness with different 
random maps, allowing one to consider more than 100 
unique arm cavity pairs. From these data, the mean and 
standard deviation of interferometer contrast defect were 
found as a function of mirror aberration RMS. Results 
for all three beams are shown in Fig. 9. The simulated 
contrast defect for Gaussian beam (TEMOO) is consistent 
with the measured value of LIGO, which is around 10~^ 
[25] and low enough for the effective detection. 
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Figure 9: Interferometer contrast defect as a function of test 
mass surface roughness (all 4 mirrors are perturbed with ran- 
dom phase map at the same level of RMS). Solid markers re- 
port mean values of numerical results with the corresponding 
shaded regions illustrating one standard deviation (see IV B), 
which is roughly four times higher than the trace in Fig. 7 
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Figure 10: Analytic calculation with different conditions for 
reducing the contrast defect: solid curve is the original con- 
trast defect; dashed line has corrective rings added to the 
phase map; dotted curve is with detuned injecting laser fre- 
quency. 



A. Better polishing 

The most direct approach is to reduce the mirror fig- 
ure error. However, reaching appropriate levels of surface 
roughness is beyond the capabilities of current technol- 
ogy. We estimate that in order to achieve reasonable 
performance, LGs^a modes require mirrors with an RMS 
roughness roughly one order of magnitude smaller than 
is currently achievable (assuming the mirror coatings in- 
troduce no additional roughness, i.e. perfectly smooth, 
uniform coatings). 

In the remainder of this section, we thus consider more 
unconventional means of reducing the contrast defect. 
Fig. 10 shows the results of each investigation. 



B. Arm cavity detuning 



Our numerical work confirms the result from pertur- 
bation theory; LGa^s interferometers are more sensitive 
to mirror surface roughness than those supporting a fun- 
damental Gaussian mode. We further show that LGs^a 
beams are also outperformed in this respect by nearly- 
concentric mesa beams, indicating that this sensitivity 
arises due to the properties of the Laguerre-Gauss mode 
itself and is not an inevitable handicap for all beams ca- 
pable of mitigating mirror thermal noise. 

V. CONTRAST DEFECT IMPROVEMENT 

Here we examine several methods of reducing the con- 
trast defect. 



Equation (12) shows that the output field varies with 
the frequency of light injected into a cavity, or equiva- 
lently with cavity length. This motivated us to study 
the variation in contrast defect as a function of arm cav- 
ity detuning. 

It was found that detuning was effective in modifying 
contrast defect. Unfortunately, the large detuning nec- 
essary to recover good contrast had the effect of simul- 
taneously reducing the optical power circulating in the 
cavity. 

With a detuning of 100 Hz (approximately two cavity 
linewidths), it was possible to recover acceptable contrast 
(dotted line. Fig. 10). However, the same detuning causes 
the circulating power to drop by ^60%. We therefore do 
not consider this approach to be viable in gravitational- 
wave interferometers where the injected frequency is usu- 
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ally tuned to maximize the optical power circulating in 
the coupled- cavity system. This technique can, however, 
be considered for experiments where the thermal noise 
needs to be reduced, but the cavity's stored power is not 
of concern. 



C. Mirror Corrections 



reducing the stored power and thus the interferometer's 
phase sensitivity. We hence conclude that this approach 
is not promising. 

At Caltech, Yamamoto has studied a similar approach 
whereby the mirror reflectivity is set to zero at the nodes 
of the LG3,3 mode [26]. This technique was found to be 
similarly unsuitable for application to gravitational wave 
interferometers . 



The increase in contrast defect observed when using 
LGs^s modes in the presence of realistic surface roughness 
results from the presence of multiple pseudo-degenerate 
higher order modes. Here we attempt to see if this effect 
can be mitigated by depositing corrective structures on 
the mirror's surface. 

By introducing material at the nodes of the desired 
LGa^a mode it was hoped that the unwanted modes from 
the same sub-space could be suppressed. 

As a concrete example, two Gaussian rings were added 
to the random phase map at nodes 1 and 3 of the LGa^s 
mode. Each ring was of the form: 



A 

20^ 



(r- 



_g 2(R/100)^ 



(20) 



where Vp indicates the position of the different nodes. 
The frequency split is plotted in Fig. 11. Compared 
to Fig. 5, the frequency splits under this condition are 
much larger, therefore the other degenerate modes will 
be harder to excite. 
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Figure 11: Frequency shift in LG modes introduced by adding 
rings to the reference phase map. 

The analytically computed contrast defect for this case 
is plotted as the dashed line in Fig. 10. Although the de- 
fect is improved for values of surface roughness similar to 
that which is currently achievable, we find that, in order 
to significantly break the modal degeneracy, the height of 
the rings must be increased to such a degree that the in- 
duced scatter becomes unacceptably high (~ 500 ppm). 



D. Mode Healing 

Previous work [27] has shown that the presence of a 
signal recycling cavity can substantially reduce the con- 
trast defect in the case where the resonant mode in the 
interferometer is TEMo,o- The higher order transverse 
modes are not resonant in the signal recycling cavity and 
are therefore suppressed. In the LGa^s case, however, the 
signal cavity is resonant for the LGs^s mode as well as 
all of the modes which are in the degenerate sub-space. 
Therefore, our expectation is that there would not be a 
mode healing effect when using any higher order mode 
which can be split in this way. In the case where the 
signal recycling cavity is detuned to amplify the gravita- 
tional wave response at a particular frequency, the situ- 
ation could be significantly more complicated due to the 
frequency splitting shown in Fig. 5. To quantitatively 
explore the effect of the compound cavity on degenerate 
modes, further analytic and numerical work is required. 



VI. CONCLUSIONS 

In this paper, we use numerical analysis as well as per- 
turbation theory to analyze the modes of a Fabry-Perot 
cavity resonating a LGs^s beam. We prove that with re- 
alistic mirror figure errors, the real output mode of the 
cavity will change significantly, resulting in an unaccept- 
able increase of the contrast defect. 

We also investigate unconventional corrective tech- 
niques to reduce the contrast defect. While they turn 
out to be unsuitable for quantum shot noise limited inter- 
ferometers, they may have some utility for other classes 
of cavities. For LGa^a modes to function effectively, we 
estimate that surface figure errors must be reduced to 
the order of 10~^ nm RMS to achieve the required con- 
trast defect of LIGO('^ 10~^). Such precise polishing and 
coating uniformity will likely not be available for several 
years. Using high-order Laguerre- Gauss modes in stan- 
dard spherical mirror cavities appears to be a poor choice 
in light of current technologies. 

Numerical simulations using mesa and normal Gaus- 
sian beams show these beams are not so sensitive to figure 
errors. Future effort will be directed toward the construc- 
tion of a new family of optical modes which can reduce 
the thermal noise impact while simultaneously being ro- 
bust against mirror imperfections. 
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Appendix A: Contrast Defect 

Here we show that the contrast defects defined in equa- 
tions (15) and (14) are equivalent in the limit of small 
perturbations. We denote the input and output field for 
the two cavities (X and Y) as |in)x, |out)x, |out)y . 

In our analytical calculation, we assume the input and 
output fields of each cavity are normalized, because we ig- 
nore transmissivity of the ETM, and the diffraction loss, 
so that we have 

(in|in)x = (in|in)y = (out|out)x = (out|out)y = 1. 

(Al) 

We can then write 

Px = (in|in)x = 1, Py = (in|in)y = 1 (A2) 

and the power at the anti-symmetric port can be writ- 
ten as 

Pas = |||out)x-|out)y|p 

= 2 — (out|out)xy — (out|out)yx- (A3) 

With the definition of contrast defect in Eq. (15), we have 

ex = 1- |(in|out)x|,ey = 1 - |(in|out)y| (A4) 

so that we can write the output field as 

|out)x = (l-ex)e*^^|in)x + |^x) 

|out)y = (1 - ey )e^^- |in)Y + \Sy) (A5) 

where 

{Sx\in)x = 0, {Sx\Sx) = ex 
((5y|in)y = 0, {Sy\Sy) = ey 

((^x|in)y = (^y|in)x=0 (A6) 

In both the analytical and numerical calculations, we 
assume one of the two interferometer cavities is perfect 
and the other is with mirror figure errors, so here we 
can write ey = 1, |(5y) = 0, therefore, to the first order 
approximation {ex <C 1), we get obtain 



which shows that the contrast defect defined in Eq. (14) 
is the same as the analytical definition in Eq. (15). 

When we consider two cavities both with imperfec- 
tions, if they are statistically independent, we can write 
{Sx\Sy) = 0, so that the average value of the contrast 
defect of the system is 



Pa5 = 2(ex + ey). 



(A8) 



We now show that Eq. (15) can be approximately writ- 
ten as given in Eq. (16). From Eq. (13), we have 



(in I out) = 



^ + ^("-"^)HnO(n^|in) (A9) 



^ J - i{uj - UJn') 

when uj — ujn' ^ 7, Eq. (A9) can be expanded as 



E 



1 



2i{uJ — UJn') 2{uj — UOn^Y 



(in|n')(n'|in) 



7 T 

(AlO) 

Note that when the mode frequency is shifted ujn' , the 
optical power in the cavity is maximized, so that the 
linear term vanishes, therefore that the modulated fre- 
quency of the beam uo can be given by: 



cj= ^(33|^/i|33) 



(All) 



The contrast defect defined in Eq. (15) is 



-{in|n ){n |in) 



Pas = 2ex 



(A7) 



^{in\n')ojl,{n'H-u;^'^ (A12) 



After perturbation, the frequency split is the eigenvalue 
of the matrix given in Eq. (9) and the eigenvector of the 
matrix is the real eigenmode of the cavity, thus we can 
write 

^(in|n')c^^,(n'|in) 

= E Hn'){n'\Sh\m'){m'\dh\l'){l'\m) 

= (^^^ ^(in|J/i|mO(m'|(5/i|in) (A13) 



where 7 = cT/4L. Then, in this limit, when the in- 
jected field is 1 33), the contrast defect can be written as 
in Eq. (16). 
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